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1. Introduction 


The standard model of the early Universe predicts that free electrons and protons 
combined to form hydrogen atoms at a redshift of order 1000. After this “recombination” 
epoch, almost, all matter is predicted to be in the form of neutral hydrogen. The Gunn- 
Peterson [1] test is a way of examining this prediction. In particular, the absence of a dip 
in the spectra of distant quasars on the blue side of the Lyman alpha line is convincing 
evidence that there is very little neutral hydrogen in the intergalactic medium. To get a 
feel for the conflict between the prediction and the observations, consider the limit imposed 
by Steidel and Sargent [2] on the intergalactic hydrogen density, rifj : 

Tin < 8.4 x 10 -12 /i cm -3 (1.1) 

at a redshift z = 2.64. Here h is the Hubble constant today in units of 100 km sec -1 
Mpc -1 and lies between .4 and 1. The standard model prediction [3] for the hydrogen 
density at this redshift is roughly 7 x 10 -6 cm -3 , violating the Steidel- Sargent limit by six 
orders of magnitude. 

It appears then that there is a diffuse spectrum of ionizing radiation permeating the 
Universe [4]. What is the source of this radiation? Some years ago Sciama [5] proposed that 
neutrinos may be both the dark matter and the source of ionizing radiation Neutrinos are 
particularly attractive candidates to serve this dual role because cosmology predicts that 
there are many background neutrinos; they exist in numbers comparable to rh" background 
photons. If the mass of one of the species of neutrinos is in the range 25 — 100 eV, then 
they can be the dark matter [6]. Further, if the lifetime for the process 

tUeavy — * > t'Light + 7 (1-2) 

is short enough - but longer than the age of the Universe, so most of the heavy neutrinos 
are still around to be the dark matter - then the photons coming from neutrino decays 
can fully ionize the intergalactic medium [7]. 

In this paper, we will analyze the Decaying Dark Matter (DDM) scenario and see what 
predictions it makes. In particular, for any given values of the neutrino mass and lifetime, 
we will calculate the photon spectrum today as a result of all past neutrino decays, and 
we will find the predicted value of the neutral hydrogen density. To do this, in Section 3 
we derive the Boltzmann equations which govern the photon spectrum and the ionization 
ratio. It will be shown that the ionization ratio depends on the electron temperature, so 
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we will derive the Boltzmann equation for this quantity as well. Our Boltzmann treatment 
complements and - to some extent - extends previous work on the subject [81 [9], 

The possibility of radiatively decaying particles has been raised in many different 
cosmological contexts [10] [11] [12]. In all of these cases, it is useful to have an expression 
for the photon spectrum as a function of redshift. In Section 4, we solve the Boltzmann 
equation for the photon spectrum in terms of the ionization ratio. Our discussion is kept 
general so it can be applied to different decaying particle scenarios and also to other 
proposed explanations of the Gunn- Peterson tests. 

In section 5, we present the results of a numerical integration of the Boltzmann equa- 
tions. There are two requirements which together bound the neutrino lifetime. The lifetime 
must be short enough so that almost all hydrogen atoms axe ionized. If the lifetime is too 
short, though, we would have observed some of the photons from the neutrino decays in 
the diffuse UV background. These complementary bounds constrain the neutrino lifetime 
to be of order 10 23 - 10 24 seconds. Besides the quantitative results, our most striking 
prediction is that the photon spectrum today will be sharply cut off at long wavelengths. 

First, though, for the purpose of orientation, several general comments are presented 
in the next section. 


2. General Considerations 

With very hi tie work, one can get a feel for the range of neutrino masses mu lifetimes 
necessary for the DDM scenario to work. First, we note that cosmological neutrinos with 
a mass in the desired range are very non-relativistic today. Therefore, the energy of the 
photon emitted in a decay is E v f 2 = m„/2 [as long as the emitted neutrino is much less 
massive than the decaying one]. This must be greater than the energy required to ionize 

a hydrogen atom, e 0 = 13.6 eV, so the first constraint is on the mass of the decaying 
neutrino: 

> 27.2 eV. (2.1) 

A second constraint emerges when we note that if the energy density of neutrinos is 

enough to set pj ^critical = = 1, then there is a relationship between the Hubble constant 

and the neutrino mass. This follows since 
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If we now require % + % a ryon = 1, then we can express h as a function of m„: 


h(m v ) 


' m v 
.91.5 eV 


1 1/2 


+ .0125 


(2.3) 


Here - and throughout - we have set f^b&ryon^ 2 = .0125, the favored value emerging 
from cosmic nucleosynthesis calculations [3]. We note that since m„ must be greater than 
2eo = 27.2 eV for the emitted photons to be able to ionize hydrogen, h must be greater 
than 0.55. Since we are assuming a matter-dominated Universe with ft = 1, h must be 
less than 0.65 in order for the age of the Universe to be greater than 10 10 years [13], This 
places an upper bound onm„: 

ra* < 37.5 eV. (2.4) 

A final constraint comes from requiring that all the hydrogen atoms in the intergalac- 
tic medium have been ionized. If every photon produced ionizes a hydrogen atom, then - if 
decaying neutrinos are to be solely responsible for the ionization - the number of neutrinos 
which have decayed by now must be greater than the number of hydrogen atoms. But the 
number density of neutrinos that have decayed by now is just %/r)n„%), where t 0 is 
the present age of the Universe; r is the neutrino lifetime; and n„%) is the cosmological 
number density of one species of neutrinos, roughly a third the number density of back- 
ground photons. As a first estimate, let us assume that the number density of hydrogen 
[ionized or not] in the intergalactic medium is of order the diffuse baryon number density, 
4 x 10 10 times the number density of background photons. By requiring the number 
density of decayed neutrinos to be greater than the hydrogen number density, we arrive at 
the third constraint, this one on the neutrino lifetime: 


t 0 <t < 10 % (2.5) 

where the lower limit follows from the requirement that most of the heavy neutrinos have 
not yet decayed. We will see in the next few sections that the lower limit is much larger 
than t 0 due to ultraviolet constraints [14] and the upper limit is much smaller than 10% 
due to the effects of recombination. The favored regime which emerges - 10 23 - 10 24 
seconds - is extremely interesting. The standard model of particle physics predicts [15] a 
lifetime much longer than this for a 30 eV neutrino, so convincing evidence of this scenario 
would require drastic modifications of the standard model [16] [17]. 

How good is the assumption that the number density of diffuse hydrogen today is 
of order the cosmological baryon number density? The answer to this question depends 
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upon how much and when baxyons clumped together. One might imagine that all of 
the baxyons clumped together (into pop III objects) shortly after recombination. Such a 
scenario would explain the null Gunn-Peterson result but is rather unlikely as a model 
of structure formation. It is more likely that in the era in which quasars formed (the 
epoch the Gunn-Peterson test probes) virtually all the baryonic matter is still diffuse. 
Subsequently more clumping, though how much is unknown, took place. In order to allow 
for this uncertainty we shall introduce a clumping factor / such that 

_ diffuse _ r cosmological fn c\ 

n B = J n B (2.6) 

where / < 1. It is fairly safe to assume that for the oldest quasaxs / is of order unity, 
while today it could be up to three orders of magnitude lower [18]. 


3. Evolution Equations 

In this section we derive the Boltzmann equations which govern the photon spectrum 
and the ionization ratio [19]. In a flat FRW Universe, the distribution function f(k,t) of 
a paxticle species xp evolves [20] according to the following equation: 

{f (3.1) 


where R(t) is the scale factor of the Universe. The term on the left hand side is the 
generalized Liouville operator acting upon /^, which takes into account the cosmological 
expansion. The term on the right side, represents all the physical interactions involving 
the particle. For example, for scatterings xp + a «-► b + c, the collision term is 

C (k t) = — / [ <* 3 Pfe f <PPc 

* ’ 2 E*(k)J (2tt) 3 2 E a (p a )J (2tt) 3 2 E b (p b ) J (2n^2E c (p c ) 

X ( 2 ? T) 4 S 4 (k+p a - Pb - Pc)[\M\ 2 n)+a _ b+c U(k)f a {p a ) ( 3 - 2 ) 

- \M\ 2 b+c _ a+ll MPb)fc(Pc)]- 


Here, and throughout, E t {q) = \/q 2 + mj for any particle i. The matrix elements squared 
include summations over the spin states of particles a, b and c but not of xp. Note that we 
have not included Pauli (Bose) suppression (enhancement) factors, as they are negligible. 
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We shall first consider the distribution function of the photons, f y (k, t ). The aim of 
this paper is to elucidate the most important physics of the reionization epoch, so accord- 
ingly only the dominant terms will be included on the right hand side of the Boltzmann 
equation. For the photons this includes the decay process 


v H{q) -*• t'Liph(k)] 


(3-3) 


the heavy neutrino decaying into a light neutrino [assumed massless] and a photon. In 
analogy to Eq. (3.2), we have 

r _ 1 / d 3 p f d 3 q 

decay " YkJ mJ ( W2E,(qy - P - k) (3.4) 

where F p (q, t ) is the distribution function of the massive neutrinos and |M| 2 is the square 
of the amplitude for radiative decay, summed over massive and light neutrino spin states. 

The inverse process is unimportant at the temperatures of interest, far below the heavy 
neutrino mass. 

We know that F v {q) constrains q to be of order the temperature; since we are interested 
in temperatures much less than m„, we have 


<*[£(?) ~ W - k\-k]& ~S(k - m v j 2). 


(3.5) 


Using this simplification and the fact that 


[ q 

2 J ( 27 t)3 F ^) = (3.6) 


where n„(t) is the number density of the heavy n and [hence the factor of 2]. we find 


r _ n„(t) 7T 2 c/t 
'-decay — 777 — d(k — &o). 


Tkl 


(3. 


Here fc 0 is the energy of the decay-produced photon. 


h = — . 


(3.8) 


and r, the lifetime of the decaying neutrino, is given by r = lQnmJ\M\ 2 . 
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Another important process affecting the photon distribution is photoionization of 
hydrogen - 7 (k)H(k') — ► e(p)P(p') - and the inverse process, recombination - e(p)P(p') —* 
'y(k)H(k'). These lead to a collision term 

r. 1 ( l i3 f f 

2k J (2n) 3 2Ep(p') J (2ir) 3 2 E.(p) J (2w) 3 2E H (k') (3.9) 

x |M|5 / ,_,p(2i) , i 4 (p + p' -k- k')[f.,(k)gH(k') - g,_(p)gp(p’)l 

where gH(k'),g e ,gp are the distribution functions of hydrogen atoms, electrons and pro- 
tons, and M is the amplitude for the process in which a hydrogen atom in the ground 
state gets ionized. There is no contribution to the photon spectrum due to ionization from 
excited states of hydrogen, since at these low temperatures, all excited states that might 
be produced immediately drop down to the ground state: only the n = 1 state is occupied. 
[Appendix A treats this question quantitatively.] All the integrals in Eq. (3.9) can be 
carried out assuming thermal equilibrium amongst free electrons, protons, and hydrogen 
atoms. Since the rate for elastic scattering - the process which bring the species into 
thermal equilibrium - is much greater than the expansion rate, thermal equilibrium is a 
good assumption. Consider first the term in Eq. (3.9) with / 7 < 7 h . The integrals over p 
and p' can be rewritten in terms of the cross section for photoionization, and the integral 
over k' can be rewritten in terms of tin yielding, 


^ionization — ^'H^Hy~*ePff{k,t')Q{k Co)- (3.10) 

Now consider the term in Eq. (3.9) with g t gp. Due to thermal equilibrium cicb of these 
distribution functions axe Maxwellian with a common temperature T e : 

“ T ' Mp,) “ - (ilif 

where M is the proton mass and the factor of 2 in the denominator for g e accounts for the 
two spin states of electrons [nucleax spin states play no role, since those in gp cancel those 
in gtf]. Note that the proton distribution function, gp , is proportional to np, which we 
have set equal to n e since there are an equal number of free electrons and protons. [Any 
difference between the two is due to the small fraction of helium atoms, which we neglect.] 
These expressions for the distribution functions together with energy conservation can be 
used to do the integrals in Eq. (3.9). Specifically we write 


9e(p)gp{p') = y 


2tt 

m,T, 


3/2 ( 27T 
\MT t 


n: 


t IT 


m e T e 


3/2 


( 2n V /2 ggjdii (ML 

\MTj n H v 2 1 


,-[k-t 0 + k' 2 /2M]/T e 


3/2 


(3.12) 
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where the first equality follows from energy conservation and the second from the fact that 
gH is also assumed Maxwellian. We can now proceed as before, integrating over p.p', and 
k' to find 

^recombination = tf 'r-ePe _[ *“ <oJ/Te . (3.13) 

Physically we see that at low temperatures the photons produced in recombination will 
have energies very close to e 0 , as expected. 

While there are other processes which in principle should be included in the Boltzmann 
equation, such as Compton scattering and Bremsstrahlung [21] , since the rates of these 
processes are much slower, they can be neglected. Thus the photon equation is 


/ _ R i.JL\f(u.k\ n„(t)Tr 2 

(dt R k dk) f ^ k, ^~ Vk 2 6 ( k ~ k °) 

( 27 r \ 3/2 

+ \^T e J n ^-P e_t *- <o]/Te - n H (TH^ePMk,t)e(k - eo). 

(3.14) 

We must now derive an equation for the neutral hydrogen density. Following the 
principles above, we can write the non-integrated equation as 


. 1 / f d y r d 3 k 4 

2 E H (k’) J (2tt) 3 2 E e (p) J (2tt) 3 2 E P (p') J (2tt) 3 2A: ^ 27r ^ 6 (P + P 1 - k- 

X [l M\\ p _„^'(p)g P (p') - \M\* Hy _ tP Mk)g„(k')} 


k’) 


(3.15) 


where we have explicitly distinguished the two amplitudes for reasons to be made clear 

shortly Summing over hvdrogen spin states and integrating over hydrogen momenta we 
find 

R ~ 3 J t ( n H R3 ) = 2n P J -^~9e{p)(a eP ^ Hl v) 

f dH- (3 - 16 > 

~ 2n H J 

The recombination cross section a eP ^ H ^ must include capture into all possible states of 
hydrogen, since at low electron temperatures, direct recombination to excited states gives 
an important contribution to the recombination rate. Since the lifetime of excited states 
(even metastables) is much shorter than the time scale of ionization, a H ^ eP need only 
include ionization from the ground state. 



To simplify the equation further we define the neutral hydrogen ratio 


_ n H 


n B 


(3.17) 


where ns is the number density of baryons in the intergalactic medium. We are neglecting 
the small fraction of helium atoms, so n e =-n p = ns( 1 — r). Plugging in the relevant cross 
sections for ionization and recombination [22] and performing the integrations, we have 

^ = n s (cru)fl(eo/Te) 1/2 (0.43 4- l/21n ^-)(1 - r) 2 - n 7 cr/r (3.18) 


where 



(3.19) 


Here <jj = 6.2 x 10 -18 cm 2 is the ionization cross section at threshhold, and we have taken 
t TH-y—tP to fall off as k~ 3 . The recombination cross section is (<tv)r = 5.2 x 10 -14 cm 3 
sec -1 . 


Since the recombination rate depends upon the electron temperature, we must derive 
an equation for it. We begin with the Boltzmann equation for the electron distribution 
g e (p,t). There nr** several relevant interactions which must be included on tlv right hand 
side of the Boltzmann equation. First, there are elastic processes: Coulomb scattering 
with other free electron; Coulomb scattering with free protons; and elastic scattering off 
hydrogen atoms. We also must include a term representing recombination and its inverse, 
photoionization. Finally, a free electron can scatter off the background photons and lose 
energy, so we must include a term for Compton scattering. Thus the equation for the 
electron distribution function is 


d R d 


fa 


~^P- 


■ loA n. t ) 


We now outline a simple recipe for dealing with the elastic collision terms, (i) Multiply 
this equation by the kinetic energy p 2 /2m e and integrate over phase space: number of 
spin states xd 3 p/(2Tc) 3 . (ii) Do the same for all particles with which the electron interacts 
elastically, (iii) Add all these equations. Then the integrals over the elastic terms all vanish 
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due to energy conservation. As an example, consider electron -proton elastic scattering. 
The sum of these - after the relevant multiplications and integrations is - 


/ |^/2m,)C + / 0ji(p' 2 /2M)CZU Uc 


- f f <Pp' f d 3 q f d 3 q' 

J (2tt) 3 2 E e (p)J (2ny2E P (p')J (2iry2E t (q) J (2tt) 3 2 E P {q') ( 3>2 1 ) 

x (27r) 4 6 i (p + p' - q - q')^2 l-^l 2 


spins 


x [p 2 /2m t + p' 2 /2M][g e (q)gp(q') - g e (p)gp(p')}. 
But, by energy conservation, the last line here is 


(? 2 /2m e + q' 2 /2M)g e (q)g P (q') - (p 2 /2m t + p' 2 /2M)g e (p)g P (p'). (3.22) 


This changes sign under the interchange of dummy variables (q +-► q';p <-► p') while the 
integration measure is invariant under such a change. Therefore the sum in Eq. (3.21) 
vanishes. 

The above prescription throws away some of the information contained in the non- 
integrated equation. However, we can compensate for this by assuming that the distribu- 
tions are Maxwellian. Then the integrated equation determines the one parameter in the 
Maxwellian distribution: the temperature T e . Specifically, it is straightforward to show 
that the left hand side of the sum of the integrated Boltzmann equations is 

+ < 3 - 23 > 

We can write the densities in terms of r and np and then use the fact that npR 3 is constant 
to rewrite this as 


\n B R- 2 j t {R 2 T t (1-r)} = ln B {2 




(3.24) 


If the right hand side of the Boltzmann equation were zero - that is, in the absence of 
any interactions - and if r were 0 then we find that T e falls as a well-known result. 

Another good check on Eq. (3.24) is to ignore expansion and interactions. Then as the 
neutral hydrogen fraction increases (r > 0), T e increases as well. This reflects the fact that 
the temperature is shared by the 3 species: electrons, protons, and hydrogen atoms. If one 
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hydrogen atom is added, both an electron and a proton must be subtracted, so the total 
number of degrees of freedom has dropped by one. This leads to a rise in temperature. 

It remains T o find the right hand side of the Boltzmann equation; that is. to integrate 
the remaining collision terms in Eq. (3.20) and its proton and hydrogen counterparts. 
First consider the recombination/ionization terms. Schematically they sum to 

X l M l 2 - 9e{p)gp(p')} [(p 2 /2m e ) + (p' 2 /2M) - ( k n /2M )]. (3.25) 

momenta, spins 


By energy conservation, the energies in square brackets are just equal to k - e 0 . Thus, 
when integrating over f y gH, the p and p' integrals can be done easily. When integrating 
over g e gp , we can simply set the energies equal to the electron kinetic energy, since the 
heavy atom gains very little kinetic energy in the process. Therefore, the k and k' integrals 
in this case are easily performed. The final integral over the electron momentum is just 
the thermally averaged cross section multiplied by the temperature. So these two terms 
give: 


/ d 3 k 

^^/ 7 (fc)(fc-eo)(<Tff7-ep) - n e n P T e ((7t;) fl (eo/r e ) 1/2 (0.43 + l/21n^-). (3.26) 

Finally the Compton term in Eq. (3.20) must be integrated. Note that Compton scattering 
off the background photons is not important for protons [and obviously not for neutral hy- 
drogen] since the cross section is small than for electrons by a factor of ( m e /M ) 2 . Therefore 
we need only integrate over Ccompton of electrons. The result is well- known |23] 


2 J (^(p ! /2m < )Ccon,p.o» = Zn,iT 2 ~ — 


(3.27) 


where T is the temperature of the background radiation and t c is the Compton cooling 
time 


TV = 


3m e 

4p*ir ^Thompson 


= r c °(T 0 /T) 4 


(3.28) 


where the Compton cooling time today (when T = T 0 ) is r° = 7.1 x 10 19 sec. We can 
collect equations (3. 24), (3. 26), and (3.27) to write the evolution equation for the electron 
temperature as 



I* 1 - + 4 t ‘ - ftt 7 "'} = 2 7 

+ |(1 - r) - re (r/T 0 ) 4 - (1 — r ) 2 np((Tt;)p(eo/Te) , ^ 2 (0.43 4- l/21n ~ )T e . 

T c T e 

(3.29) 
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4. Photon Spectrum 


The set of equations (3.18), (3.29), and (3.14)- can be numerically integrated to 
find the photon spectrum and the ionization fraction at all times. We will display results 
of this numerical work in the next section. Here we focus on the equation for the photon 
spectrum - Eq. (3.14). For a variety of reasons it is useful to have an analytic expression 
for the occupation number, / 7 (t/,z). We first derive such an expression and then apply it 
to three problems that arise when dealing with a non-equilibrium distribution of photons 
in the early Universe: (i) What is the remnant spectrum of photons today? (ii) How many 
ionizing photons are present at a given redshift z? and (iii) Why and when is recombination 
to the ground state of hydrogen suppressed? This last point - initially made by Peebles 
will be shown to emerge in a straightforward manner from the analytic expression for 

My,z)- 


It is useful first to rewrite the photon equation in terms of dimensionless variables. 
Let 


y = 


A 

k o 


1+z = 


R(to) 

m 



(4.1) 


where the last equality follows from our assumption that the Universe is flat and matter 
dominated. By way of orientation in these new variables, we note that all photons are 
produced with energy k = k Q and lose energy due to the redshift, so we are interested only 
in y < 1- TJi e mam redshifts z of interest are determined by the neutrino lifetime; roughly, 
we will be studying the epoch z < 300. Using the facts that (i) the Universe is matter 
dominated and therefore f = J?o(l + *) 3/2 (H 0 is the Hubble constant today) and (ii) the 
baryon and neutrino number densities increase as (1 + z) 3 , the photon equation becomes 


3?r 2 ton„(to ) 


2 rk 3 


-S(y - 1 ) (1 + ~) 3/ 2 


(4.2) 


+ A 


2tt(1 + z) 3 \ 3/2 


m e T e 


n B (t 0 )(l - r ) 2 e- ( * oj '- e ° ) / ;r ' 0 [y - €o /jfc 0 ] 


where 


_ n B {t o)< 7 / 

Ho 


(4.3) 


The first term on the right in Eq. (4.2) is the source term due to decays and the second 
the source term due to recombination. The effect of the ionization sink term is governed 
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by the dimensionless number A, which is essentially the ratio of the ionization rate to the 
Hubble rate. We can get a feel for the magnitude of A by taking riB(to) = 2 x 10 -7 cm -3 
- the favored value from cosmic nucleosynthesis. Then the ionization rate today [if all 
baryons were in the form of neutral hydrogen] is n B (to)cri = 3 x 10 -14 sec -1 This is much 
larger than the expansion rate today, so A is a very large number. 

For the time being we will call the right hand side of Eq. (4.2) S(y, z) so that the 
final form of the solution may be used for any sources of interest. We note, though, that S 
typically has [at least] two pieces: the source of ionizing photons - in this case the decaying 
neutrinos - and the photons produced in the process of recombination. To solve this partial 
differential equation, we first turn it into an ordinary differential equation by introducing 
a new variable 


p = — In y — — In 


1 + z 

l + z c 


(4.4) 


The variable p runs along a curve in the (y, z) plane; different curves are labelled by 
different values of z c . Along any one curve, 


d_ 

dp 


-(1 + *) 


dz 



so that in terms of p , Eq. (4.2) becomes 


(4.5) 


#vCg) 

dp 


+ MM/hf -I (1 + *‘>^1 IS + 3 A1^£L - «/*»> 


— S [e (1 + z c )e pl \ 


(e-p)3 




(4.6) 


To be explicit, the right hand side here means: take S(y, z) and evaluate it with y set equal 
to c p and 1 + 2 set equal to (1 + z c )e p . For a given z c , the solution to this ordinary 
differential equation is elementary: 


f-Ap) 


-s: 


dp's 


e p , (1 + z c )e 


-p 


f fP A(e 0 /fco) 3 r 

v PYn < 1 A r\ U 

(1 + z c )e p " 

(1 + 2 c )e p 

3/2 

0| 

O 

0 
yj 

1 

1 

J P . " e-3p” J 


Here we have dropped the boundary term, f^ip = 0), since there are no photons with 
momentum greater than &o, i- e - with p negative. This solution in terms of p can be 
converted back into a solution / 7 (y,z) by eliminating p [via p = -lny] and 2 C [via 1 + 


12 


Zc ~ (I "b z )!v\- It is also convenient to introduce new dummy variables: y' = e~ p ' and 
I + z> = (1 + Zc)e~ p = . These steps lead to 


A(V.*)= f ^-Sly',(l + ,) y -) 

Jy y y 

( , /i + ?\ 3 / ,(1+r )»'/» 

x exp l - A(e 0 Ao) ( ) / 

I \ V / imaiJl + tjfl+ilji 


dz' 

] (l+z') 5/2 


(4.8) 


r(*') l. 


We can go no further without an explicit form for the source terms, but we note that A 
is a very large number - of order 10 3 - so f y (y, z) is negligible unless other factors in the 
argument of the exponential are very small. That is, ionization - as represented here by 
the argument of the exponential - plays an important role in suppressing the number of 
photons coming from any source. 

In our case there are two source terms; we focus first on the dominant one: the photons 
coming from decaying neutrinos. The contribution to / 7 (y,z) from this term is 


■>(»**•'■) 


3/2 


I , / 1 + z\ 3 + 

x exp - A(e 0 Ao) 3 ( ) / 

\ y / Jmaz[\ + z,(l + *) 

37r 2 t 0 n v (<o) r vV , , 

2 7^ (1 + 2) 0 ( 1 -^) 


dz 1 


(1 + z') 5 / 2 


r(z') 


(4.9) 


x exp < -A(eoAo) 


3 , i + z\ 3 r( 1+z »* 

J imai[l + j,(l+i)j 


dz' 


(1 + z') 5/2 


r(z') 


where the superscript on / 7 refers to the fact that this is the contribution from the first 
source term. Eq. (4.9) tells us that the photon spectrum at a given energy y and epoch z 
depends on an integration over the amount of neutral hydrogen at redshifts larger than z. 
This is reasonable: In the absence of ionization, high energy photons red-shift down and 
add to the number of low energy photons at later epochs. If these high energy photons 
are lost through ionization, the low energy number at later epochs is correspondingly 
suppressed. The integral over r(z) quantifies this suppression. In short, if we are given 
r(z) we can do the integral and determine the photon spectrum. In practice things are 
not so simple: r(z) depends on f y . Nonetheless, we will see that in several important 
cases, one can make good approximations which considerably simplify the way in which 
/-. depends on r. 
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We can use Eq. (4.9) to estimate the photon spectrum today when 2 = 0. Since 
2 = 0, the integral in the argument of the exponential in Eq. (4.9) ranges from z' = 0 to 
1 + z' = 1/y — ko/k. As a first approximation to the effect of ionization, let us assume 
that the Universe is spontaneously ionized at 2 = z*; i.e. r = 0 for all 2 < 2 ,. Then the 
argument of the exponential vanishes as long as 1 + Z{ > ko/k. For energies k less them 
Jfc 0 /(1 + 2 ,), the argument of the exponential is huge and'/ 7 ~ 0. We can write for the 
photon intensity today 

dfiy 2 d 3 k / 7 (fc,z = 0) 

= k die dtt = (2tt) 3 dk dO, 

3/2 (4.10) 

= 3nMto/k\ _ k) _ + 2j)) 

87r r \ko / 


The first part of this expression is the standard formula for the intensity due to a radiatively 
decaying relic particle [10]. Only the last step function - reflecting the ionization process 
- is new. Of course the drop in intensity at ko/(l + Z { ) is not as dramatic as the step 
function indicates; in the next section we will show numerical results. 

Another question of interest is : At a given redshift 2 , how many ionizing photons are 
present due to the decaying neutrinos? To find this, we must calculate n 7 , as defined in 
Eq. (3.19), 




( 2 ) -• */(t< )) 


3t 0 eo(l + z) 3 / 2 


2rk 3 


«o/*o 


e 0 (l + *)\ 3 [ (1+z)/y d(l + z') 


I -^fL 


(4.11) 


(I + 2') 5 / 2 




One simple approximation which is often useful is to neglect completely the effect of 
ionization on / 7 ; this corresponds to setting the exponential in Eq. (4.11) to 1 so that 


n 


( 1 ) 


(z) 


n 


'° nized = n„(* 0 )-(l + zf> 2 (e 0 /k 0 ) 3 


(hie o) 3/2 -l 


(4.12) 


where the superscript reminds us that this approximation should be valid only in an ionized 
Universe. For if there are very few neutral hydrogen atoms, then the ionization process has 
a negligible effect on the photon distribution. Note that as k 0 approaches eo, the number 
of ionizing photons goes to zero: expansion causes the decay-produced photons to lose 
energy, pushing them below threshold. 
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The opposite situation occurs when there are many neutral atoms so that a photon 
produced in a neutrino decay immediately ionizes a neutral atom. This corresponds to 
a large damping factor in the exponential of Eq. (4.9). In such a situation there is no 
time for a photon to lose energy via the red-shift before it is absorbed in an ionization 
process. Therefore, there are almost no photons with energy much less than kq. We can 

take advantage of this fact by expanding the argument of the exponential around y = 1. 
This leads to 



3f 0 (l + z) 3 / 2 
2Ar 


A(«o/Ao) s [l— to/fco] 
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dv 
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4 A(e 0 /*o) 3 2^1n(l + z)^A(e 0 /ifco) 3 + 0(v 2 /\ 2 ) 


din r 


) 


(4.13) 

where the dummy variable i/ = A( £o Ao) 3 (l - y). The upper limit of this integral is quite 
large, effectively oo. It is clear that as long as r(z)(l + z) 3 / 2 is of order 1, the dominant 
contribution comes from v as 1. Hence the terms of order v/\ do not contribute in this 


pre-ionization regime. Only when r(z) becomes very small (= 0(A -1 (1 4- z) -3 / 2 )) do the 

higher order terms in the argument of the exponential become significant. Hence as long 
as r(z) is not too small, 


hi^Vz) 


n„(t 0 )<o 3(1 +z) 3 / 2 
r 2A 


OO 

j d v exp| — r(z)i/(l + z) 3 / 2 J. (4.14) 


The remaining integral is trivial so that the un-ionized approximation is 


*?(*) 


^un— ionized ^^(^ 0)^0 3 

u y ~ 


(4.15) 


neu- 


r 2A r(z) ’ 

Until now we have focused on the source term for photons coming from decaying 
tnnos. Now let us consider the source term due to the process ep -> Hj. The contribution 
of this process to the photon occupation number can be obtained by inserting the last term 
on the right hand side of Eq. (4.2) into our general expression for /-, (Eq. (4.8)): 

3 + y > / 




exp^ -A(eoAo ) 3 
3/2 


)L 


maztl + c.d + cJjnL.] (1 + z')*/' 2 


•-0 
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(4.16) 
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This term can be shown to be much smaller than fy l \ However, it does have one important 
implication for the recombination process. Consider the contribution of this term to n 7 , 
the number density of ionizing photons. We show in Appendix B that, if the electron 
temperature is small, 





(1 - r) 2 0.8(<7u);t 


07 


1 + 




A(l+z)»/3rT. 


(4.17) 


This contribution to n 7 must be included in the evolution equation for r (Eq. (3.18)), the 
right hand side of which is 


nB{<rv) R (e 0 /T e ) 1/2 (OA3 + l/21n^-)(l - r) 2 - n 7 <r/r 
= n B (<ri;) fl (eo/T«) 1/2 (0.43 + l/21n - --- 5 * * 8 <n )(1 - r) 2 - n^c r/r 

le 1 + A(l + z) 8 /*r T. 


(4.18) 


The new term - coming from n 7 - reduces the recombination rate. In fact, if A(1 + 
z) 3 / 2 rT e /e 0 is very large, the factor of 0.8 exactly subtracts out recombination to the 
ground state. The fact that recombination to the ground state is strongly suppressed 
was first noted by Peebles [24] in the context of the standard cosmological scenario and 
recently emphasized by Asselin etal. [9] in the context of decaying particle scenarios. 
The reason for the suppression is straightforward: Any recombination to the ground state 
produces an ionizing photon. If the ionization time is very short compared o the Hubble 
time - i.e. if A(1 r z) 3 ^ 2 rT e / €0 is very large - then this photon will quickly ionise another 
neutral atom. Thus the total ionization fraction remains unchanged. Only when r becomes 
very small does the ionization time become larger than the expansion time; at that point 
recombination to the ground state also becomes relevant. 


5. Numerical Results 

In this section, we present numerical solutions to the coupled equations -(3.18), (3.29), 

and (4.8). For definiteness, we will first present results and discuss the details of a 30 eV 
neutrino with a lifetime of 10 24 seconds. The qualiatative behavior is independent of the 
exact neutrino parameters in the parameter space of interest ( 27.2 eV < m„ < 37.5 eV 
and 10 22 sec < r < 10 26 sec), though, obviously, the exact quantitative results will differ. 

First, we shall consider the evolution of the ionization state of the Universe. In Figure 
1 we have plotted the neutral hydrogen fraction, r, versus 1 + z. The most striking feature 
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of this plot is the precipitous drop in r - four orders of magnitude - while z changes by one 
or two. This suggests that reionization was spontaneous at a z t of about 17. If we plot the 
complementary quantity, 1 — r, a different picture emerges. In Figure 2 we have plotted 
the free electron ratio, 1 — r, versus redshift. The salient characteristic of this plot is the 
gradual change in 1 — r until total ionization. Is the ionization spontaneous or not? That 
depends upon what question is being asked. If one’ is interested in how often CBR photons 
scatter off free electrons, then the important scale is the scattering rate of the photons 
which is proportional to n t ^Thompson) where n e = (1 — r)ng. The opacity — roughly the 
scattering rate divided by the Hubble rate - changes gradually, as suggested by Fig. 2. 
This is the issue analyzed and the conclusion reached by Scott, Rees, and Sciama [25]. 
On the other hand, we are interested in how reionization affects the spectrum of decay 
photons. The relevant scale is ngai , where nn — rng • The relevant plot is Fig. 1 in 
which it is seen that the opacity changes very rapidly for these photons. 

The scale which characterizes the ionization rate is n 7 < t/. In Figure 3 we have plotted 
the ratio of ionizing photons to baxyons sis a function of redshift. The two solid lines 
correspond to the two sources of ionizing photons: decaying neutrinos (Eq. (4.11)) and 
recombinations (Eq. (4.17)). This plot corroborates the claim made in the last section 
that the recombination source is always less important than the decaying particle source, 
and thus its contribution to the photon spectrum is negligible. The importance of this 
term lies in the correction it provides to the recombination rate. In the pre- ionization 
regime, when tns ratio of these photons to neutral atoms is much smaller than uni tv, every 
recombination photon immediately reionizes a neutral atom. Thus direct recombination to 
the ground state is suppressed. However, in the post-ionization regime direct recombination 
to the ground state is unsuppressed. It is reassuring that in our formalism this comes out 
naturally. 

The shape of the recombination photon density curve provides insight about the var- 
ious mechanisms driving recombination. Early on, when the Universe is nearly neutral, 
there is very little recombination, and hence n 7 ^ is very small. As the Universe becomes 
reionized recombinations become more common, so there is a rise in n 7 2> . Eventually, due 
to the Hubble expansion, the free electron-proton density turns around again, and the 
number of recombination photons correspondingly drops. The rise in n^ 2> at z ~ z x corre- 
sponds to the physical fact that recombination [to all levels] is an important process at this 
epoch. As an example of this, we can estimate what z x would have been in the absence of 
recombination. Then, once the number of photons produced by neutrinos was equal to the 


17 


number of baryons, the Universe would be ionized. That is, a simple estimate would give 
n v ti/r = ng, or Zj = /tjt) 2/3 . With the parameters under consideration, this 

estimates z, ~ 45. The fact that z, is actually closer to 20 illustrates that recombination 
delays the onset of full ionization [26]. 

We have also included in this figure the two approximations to discussed in the 
last section. The dashed curve corresponds to Eq.(4.15), 'and the dotted curve to Eq. 
(4.12). It is evident from these plots that each approximation is extremely good in a 
particular epoch. The first approximation is virtually exact up until z « z*. The reason 
this approximation is so good in the un-ionized regime is straightforward. Let’s rewrite 
Eq. (4.15) in the more transparent notation, 


^(l) ^ T M<)/ r 

fly 


(5.1) 


This says that the number of photons is approximately equal to the rate of neutrino 
decay divided by the rate of ionization. When the ratio of photons to neutral atoms 
is so small that every photon will immediately ionize a neutral atom then the Hubble 
expansion is unimportant, and the approximate equality is virtually exact. With our 
sample parameters it turns out that even at z = z,, there are about 160 neutral atoms for 
each decay photon, so it is not suprising that the approximation is so good up until this 
point. It is only when the number of neutral atoms per photon is of order uv.> , m d less, 
i.e. for z < z,, that approximation (5.1) breaks down. We see here that this r, ?,n extremely 
good approximation until r falls below 10~ 4 . If one is interested in calculating the density 
of free electrons then this approximation is sufficient. However, if one is interested in a 

calculation of the neutral atom fraction for z < Zj, then this approximation breaks down 
too soon. 

The post-ionization approximation, Eq. (4.12), becomes accurate a few redshifts 
after the un-ionized approximation becomes invalid. Recall from the discussion in the last 
section that this approximation neglects ionization totally. In Fig. 1, we see, however, that 
r continues to drop even in the post-ionization regime. The residual fraction of neutral 
atoms is not strictly frozen, but rather continues to get smaller, so there still is some 
ionization happening. Since the ratio of the ionization rate to the Hubble rate decreases 
at lower redshifts, this approximation becomes virtually exact. Nonetheless, there is a 
regime, albeit a fairly small one, in which neither approximation obtains, so numerical 
integration of the photon spectrum, eq (3.19), is critical. 
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The DDM scenario predicts a diffuse spectrum today. Armed with the ionization 
history of the Universe in the reionization epoch, we are able to calculate it, by numerically 
integrating Eq. (4.9) with z = 0. Others who have considered the DDM scenario [10] have 
ignored this effect . 

In Fig. 4 we have plotted the correct DDM spectrum (solid line). For comparison’s 
sake we have also plotted the uncorrected spectrum, i.e. the one which ignores ionization 
effects (dashed line). The two spectra are identical for photon energies greater than a 
few eV. At a low photon energy (for our sample parameters at about k = .9 eV ), the 
DDM flux drops quite abruptly to zero. This is readily understood in light of previous 
discussions. At z ss z, there is a precipitous drop in opacity. Photons today with energies 
^ < p+^7 would have been emitted while the Universe was still opaque, i.e. when z > z,. 
Essentially every photon produced in this pre-ionization regime is used up in the process 
of reionization. Photons today with energy k > were produced in a transparent 

Universe, so accordingly travel freely. In the last section this was approximated as a 
step function at k = &o/(l + Zi). The numerical result verifies that this is a reasonable 
approximation. 

Fig. 4 also shows several sets of data points. The two in the low energy regime 
[27] are representative of upper limits in the cosmic flux: no detection has been made. 
Any predicted cosmic flux must be lower than this. It is seen that the upper limit is fax 
above the predicted levels even in the absence of ionization. Therefore, while data in this 
regime may ultimately be used to detect the radiation emitted by unstable neutrinos and 
in particular the characteristic drop in intensity due to ionization by early photons, the 
present upper limit is several orders of magnitude away from this goal. 

The second set of data points in Fig. 4 lie between 6 - 10 eV [28]. They purport to be 
actual measurements of a cosmic background flux with all local contaminants subtracted 
off. These measurements are much closer to the spectrum produced in the DDM scenario 
for two reasons: first, the observed magnitude is an order of magnitude or so lower than 
the upper limit at 3 e\ and, second, the predicted flux is larger here, since it grows as 
^ ^ ■ Therefore, this regime is most likely to be of use in detecting photons from unstable 
neutrinos. The only unfortunate aspect of this is that the sharp drop in intensity due to 
ionization occurs at lower energies than this most favorable regime. 

Since the predicted spectrum is a function of (m„,T), the observed points can be 
turned into constraints on the neutrino mass and lifetime. We wfill see shortly that this set 


19 


of constraints complements the constraints coming from the Gunn-Peterson test to allow 
only a small window in parameter space. 

For completeness we include a plot of electron temperature versus redshift (Figure 5). 
It is seen that the temperature does not change very drastically. In fact the temperature 
remains low enough so that collisions] ionization can be justifiably neglected. We also note 
that the thermal evolution is insensitive to the initial temperature. Compton scattering 
insures that the electron temperature tracks the photon temperature until the decays set 
in. 

The second prediction of the DDM scenario is that there is a relic abundance of neutral 
hydrogen. When Sciama and others [5] first considered the DDM hypothesis in light of 
Gunn-Peterson, they required simply that total reionization occur by z — 4.7 since this is 
the epoch of the earliest test [29]. As discussed earlier, the approximations they made do 
not allow a precise calculation of the neutral atom abundance for z < Z{. By computing r 
all the way to z = 0, we can simply calculate 

n H {z) = r(z)n B (z ) = r(2)/riB(<o)(l + zf (5.2) 

where f is the clumping factor defined in Eq. (2.6). This is plotted in Fig. 6 for / = 1. 

The import '.•nee of this prediction is that the neutral hydrogen density ; .in be probed 
observationally via the Gunn-Peterson test. In figure 6 we have also plotted the upper 
limits on the neutral hydrogen density coming from three different Gunn-Peterson tests at 
2 = 4.7, 4.1 and 2.64 [29] [30] [2]. 

The major uncertainty in our calculation of the neutral hydrogen abundance stems 
from the uncertainty in clumping. Throughout we have assumed a homogeneous and 
isotropic Universe, but clearly once quasars, whose observation allows one to deduce the 
neutral hydrogen density in the IGM, are formed (and possibly earlier) the Universe is no 
longer homogeneous. While the data point at 1 + z = 5.7 appears to impose the least 
stringent requirement, we have included it since at higher redshifts it is safer to assume 
that / is of order unity. 

Finally, we are in a position to summarize the constraints on the mass and lifetime 
of the dark matter particle in the DDM scenario. In Figure 7, we have collected the 
constraints coming from the Gunn-Peterson tests(with / = 1) and from the UV flux data. 
We have presented the constraints coming from two G-P tests [29] [30]. The dashed line 
corresponds to the constraint coming from the z — 4.1 quasar. This appears to impose the 
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strictest constraint, leaving only region A in parameter space. In particular the smallest 
allowed mass is 28.5 eV with a 2 x 10 23 sec. lifetime. The dotted line, corresponding to 
the 2 = 4.7 Gunn-Peterson constraint, allows region B in parameter space, as well. This 
does not rule out, for example, a 27.7 eV dark matter particle (Sciama’s preferred value 
[7] ), provided the lifetime is larger than 2.5 x 10 23 seconds. We have included this less 
stringent constraint because it is the one in which our approximations - homogeneity and 
no other sources heating the electrons - are most likely to be valid. 

One feature of this plot should be explicitly noted. That is the fact that as ro„ 
approaches 27.2 eV the upper limit on the lifetime drops steeply. It turns out that in order 
to satisfy the Gunn-Peterson tests, it is necessary that some ionization take place in the 
ionized regime. As discussed in section 4, the number of ionizing photons in the ionized 
regime goes to 0 as k 0 approaches e 0 - The only way to compensate for this reduction 
in ionizing photon number is to lower the lifetime. Consequently, in order to satisfy the 

Gunn-Peterson tests, for masses close to 27.2 eV, the lifetime of the decaying particle must 
be short. 

More realistic accounting of clumping and electron heating may well modify some of 
these quantitative conclusions. We note, though, that the three sets of observations - of the 
Hubble parameter; diffuse photon spectrum; and neutral hydrogen density - complement 
each other. Measurement of the Hubble parameter (or the age of the Universe) seems to 
require the neutrino mass to be less than 37.5eV. Measurement of the diffuse photon 
spectrum constrains the neutrino lifetime to be greater than 10 23 seconds. The Gunn- 

Peterson tests complement these two by providing an upper limit on the neutrino lifetime 
and a lower limit on the neutrino mass. 


6. Conclusion 

We have quantitatively analyzed a scenariio of reionizing the Universe. The numerical 
results allow accurate comparison between theory and observation. In particular, we are 
able to derive strict constraints on the parameters of the decaying neutirno scenario - i.e. 
on the neutrino mass and lifetime - by requiring the Universe to be ionized enough to 
satisfy the Gunn-Peterson tests and by limiting the ultraviolet radiation produced to be 
under the observational limits. 

Several aspects of this work may be of use in other reionization scenarios and in more 
general work on the early Universe. 
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(i) We have derived an expression, Eq. (4.8), for the photon spectrum due to a general 
source in the presence of a background of neutral hydrogen atoms. 

(ii) A qualitative feature of the photon spectrum in such a situation is a sharp drop in the 
spectrum today. Photons with energies lower than this cut-oof would have been produced 
before the Universe was ionized. Hence, they were immediately absorbed in the ionization 
process and axe no longer present today. 

(iii) The number of neutral hydrogen atoms drops dramatically at the time of ionization. 
This is completely consistent with the point made by several groups that the number of 
free ions rises gradually. In our language r changes rapidly while 1 — r changes very slowly. 
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Appendix A. Ionization From Excited States 


Here we justify our claim that direct ionization of excited hydrogen atoms can be 
neglected. It is useful to first review the standard recombination scenario as described by 
Peebles [24]. Recombination occurs when there are essentially no thermal photons available 
to ionize ground state neutral atoms. In fact the only photons around with energy greater 
than a Rydberg come from recombination to the ground state. These photons immediately 
reionize the atoms, effectively cancelling every recombination to the ground state, such that 
neither recombination or ionization involving ground state atoms needs to be included in 
the ionization equation. The important physical processes axe direct recombination to 
excited states, ionization of excited atoms, radiative transitions between atomic levels and 
the Hubble expansion. Peebles derived the following equation to describe the ionization 
ratio, 


dr 

dt 


a f 


rn B ( 1 - r ) 2 - /3 exc r 


n erc -(10.2 eV/T) 

n H 


C 


(A.l) 
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where C is a time-dependent factor (defined explicitly in [24]) which incorporates the 
effects of the physical processes mentioned above; a exc is the excited state recombination 
coefficient, of order (o , v)fl(e 0 /T e ); /3 erc is the rate of ionization of excited atoms', and n txc 
is the number density of hydrogen atoms in excited states. The ionization rate is well- 
approximated by 173.407 where 773.4 is the number of photons with energy greater them 3.4 
eV. The upshot is that even though the the fractional population of excited states is small, 
since there is an exact cancellation between ionization and recombination of ground state 
atoms, the relevant physics concerns the excited states. 

In the reionization epoch, z < 300. all of these physical processes axe still taking place, 
so the two terms in (A.l) in principle should be included in the ionization ratio equation. 
(We note that through out this epoch, the multiplicative factor, C, is equal to unity.) 
However, in this epoch there is a source of photons energetic enough to ionize ground state 
atoms - decays. Thus to Eq. (A.l) we need to add the ionization term, —h^oin gs /nB, 
where n gs is the number density of hydrogen atoms in the ground state. 

Now let us focus on the ionization terms. Because of the exponential suppression, the 
ionization term of (A.l) - i.e. the second term on the right hand side - is clearly negligible. 
"We need to examine the new term, n^<Tjn gs /riB and determine whether ionization of 
excited atoms by decay photons needs to be included. We will be justified in neglecting 
these ionizations if the per volume rate of ionization of excited atoms is much smaller than 
the per volume rate of ionization of ground state atoms, 
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In the opaque era, z > Zi , it is easy to estimate Eq. (A.2). First of all, n 3 . 4 


~ 77 


-( 1 ) 


since the photons immediately ionize and do not get redshifted. Secondly, because there 
are virtually no resonant decay photons, the fraction of excited atoms is given by the 
expression derived by Peebles, 


n ezc _ ( A Cn tIc n 2 f \ 

n^T (r+^An ^ ) (-«> 

where A = (eo _* l)3// and A = 8.227sec _1 is the decay rate of the 2s state. The ratio 
in Eq. (A. 3) is negligibly small, always less than 10" 15 , so the condition (A.2) is easily 
satisfied. For z < s t , when photons travel freely through the universe and get redshifted, 
then 773.4 is larger than n^ 1} by up to two orders of magnitude. However, even though there 
are some resonant photons around to pump atoms from the ground state to the excited 
states, adding another term to Eq. (A. 3), the fractional population of excited states is still 
negligibly small, so criterion (A.2) is again satisfied. 
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Appendix B. Recombination Photons 

In this appendix we calculate the number density of ionizing photons coming from 
recombination. This is 

= 4 f -A 2) M 
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where we have inserted /^ 2) from Eq. (4.16). The exponential e~ ( * oS ''~ eo) / T « is very small 
except when y' is very close to e 0 Ao- But y' > y > e 0 /k 0 , so y must also be close to e 0 /fco 
and, of course, close to y' as well. We can therefore expand y' about y and expect the 
leading terms to give the largest contributions. Consider, for example, the integral in the 
argument of the first exponential: 
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Since y' — y ~ T e /c o, the higher order terms may be dropped. Keeping only the leading 
terms amouonts to setting y' = y everywhere in the integrand of Eq. (B.l) except in the 
exponentials, where only the linear terms in y' — y need be retained. Therefore,' 
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Neglecting terms of order T € /e 0 , the y- integral can be performed by setting y = e 0 /k 0 
(except in the exponential), leading to 
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The first factor in square brackets is just ng(z ) since the baryon number density scales as 
■ The second factor in square brackets is the ratio of the ground state recombination 
cross section to the ionization cross section: (<7u)o /< t/. The ground state recombination 
cross section, though, is g(eo )( (7 v)r where g is the gaunt factor, roughly 0.8 at threshhold. 
We therefore recapture Eq. (4.17). 
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Figure Captions 

1. The neutral hydrogen ratio as a function of redshift in a universe with a decaying 
neutrino of mass 30 eV and lifetime 10 24 sec. 

2. The free electron ratio as a function of redshift in a universe with a decaying 
neutrino of mass 30 eV and lifetime 10 24 sec. 

3. The density of ionizing photons as a function of redshift. The solid lines are exact 
solutions: h\ * (decay photons) is labelled (1) and fffi (recombination photons) is labelled 
(2). The dashed line corresponds to the un-ionized approximation, Eq. (4.15), and the 
dotted line to the ionized approximation, Eq. (4.12). 

4. The predicted photon spectrum due to a decaying neutrino with mass 30 eV and 
lifetime 10 sec. The solid line is the correct predicted spectrum. The dashed line ignores 
corrections due to ionization effects. The arrows mark upper limits and the crosses mark 
observations, discussed in the text. 

o. The electron temperature as a function of redshift in a universe with a decaying 
neutrino of mass 30 eV and lifetime 10 24 sec. 

6. The predicted amount of neutral hydrogen as a function of redshift in the DDM 
scenario with a neutrino mass 30 eV and lifetime 10 24 sec. The clumping factor, /, has 

been set equal ic 1. The data points are upper limits: any prediction must, be below these 
levels. 

7. The allowed values of m v and r when the clumping factor / is set to 1. The solid 
line represents the constraint coming from observations of the diffuse UV background. The 
dashed line represents the constraint coming from the * = 4.1 Gunn- Peterson test, and 
the dotted line from the 4.7 Gunn-Peterson test. Region A satisfies all three constraints. 
Region B appears to be ruled out by the z = 4.1 test, though it satisfies the z = 4.7 and 
the UV constarints. 
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